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ABSTRACT 

We present photometric observations from the Stratospheric Observatory for Infrared Astronomy (SOFIA) 
at 11.1 pm of the Type Iln supernova (SN Iln) 2010jl. The SN is undetected by SOFIA, but the upper limits 
obtained, combined with new and archival detections from Spitzer at 3.6 & 4.5 pm allow us to characterize 
the composition of the dust present. Dust in other Type Iln SNe has been shown in previous works to reside 
in a circumstellar shell of material ejected by the progenitor system in the few millenia prior to explosion. 

Our model fits show that the dust in the system shows no evidence for the strong, ubiquitous 9.7 pm feature 
from silicate dust, suggesting the presence of carbonaceous grains. The observations are best fit with 0.01- 
0.05 Mq of carbonaceous dust radiating at a temperature of ^ 550 - 620 K. The dust composition may reveal 
clues concerning the nature of the progenitor system, which remains ambiguous for this subclass. Most of the 
single star progenitor systems proposed for SNe Iln, such as luminous blue variables, red supergiants, yellow 
hypergiants, and B[e] stars, all clearly show silicate dust in their pre-SN outflows. However, this post-SN result 
is consistent with the small sample of SNe Iln with mid-IR observations, none of which show signs of emission 
from silicate dust in their IR spectra. 

Subject headings: dust, extinction — supernovae: individual(SN 2010jl) — supernovae: general 


1. INTRODUCTION 

Type Iln supernovae (SNe Iln) are named for the nar¬ 
row hydrogen emission lines produced in a dense, slowly- 
moving, pre-existing circumstellar medium (CSM). Typically 
attributed to the progenitor winds, the CSM has densities 
and velocities that suggest extrem ely high mass-los s rates, in 
som e cases up to ^ yr“' (iFox et alJl201 IL hereafter 

El 1: iKiewe et al.l[2012^ . The progenitors of SNe Iln are un¬ 
known. The implied mass-loss rates are often compared to the 
episodic dense winds observed in som e massive stars, suc h as 
Luminous Blue Variables (LBVs; e.g.. lSmith et al.ll201 IL and 
references therein). 

SNe Iln account for more than half of all known SNe with 
late-time (> 100 days) infrared (IR) emission (Ell), imply¬ 
ing the presence of warm dust. El 1 showed, via a warm mis¬ 
sion (3.6 and 4.5 pm only) Spitzer survey of 10 SNe Iln, that 
the observed emission is consistent with pre-existing CSM 
dust, heated by the optical emission generated by the for¬ 
ward shock interaction with the dense CSM. This is in con¬ 
trast to most core-collapse (CC) SNe with IR emission, where 
observations are consiste nt with a small amount of newly- 
formed ejecta dust (e.g.. iKotak et all l2005t ISugerman et al.l 
l2006t Iwilliams et alJ l2008b . Since the dust is produced in 
the outflow of the progenitor system, characterizing this dust 
may offer a clue about the pre-SN system. However, as El 1 
showed, Spitzer observations at 3.6 and 4.5 pm are sufficient 
only to establish the presence of a warm dust component, and 
not to determine the composition of the dust. 

It is well-known that massive stars make dust in their stellar 
winds, but what type of dust depends on the star’s mass and 
evolutionary state. In general, stars make either silicate or car¬ 
bon dust in their outflows, where the determining factor is the 
C/O ratio in the stellar atmosphere. The presence or absence 
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of a broad, strong emission feature at 9.7 pm that only arises 
from silicate dust provides a method by which to distinguish 
the two grain populations, provided images or spectra exist in 
the 8-13 pm range. 

1 . 1 . Potential SNe Iln Progenitors and their Dust Signatures 

The LBV scenario is a pro mising candidate for t h e pro- 
genitor systems of SNe Iln. iGal-Yam & Leonai^ (120091) 
identify a hypergiant progenitor undergoing LBV-like mass 
loss in pre-explosion images of SN 2005gl. However, it 
has several complications. Most LBVs achieve observed 
rates of only M < 10~"' Mfn y r~' in their S Doradus state 
(iHumphrevs & DavidsonllT994l) . Standard radiation and line- 
driven wind models simila rly have difficulty repr oducing 
such substantial rates (e.g., iSmith & Owockil12006^ . Clas¬ 
sical stellar evolution models suggest these massive stars 
shou ld evolve into the W olf-Rayet (WR) phase before explod¬ 
ing (iHeger et al.l l2003h . More recen tly, proposed solu tions 
include swept up Wolf-Ray et winds (!Dwarkadasll201 Ih . bi- 
nary accretion (lKa shi||2010l). and gravity-wave driven winds 
(lOuataert & Shiodel 120121) . but these all have limitations of 
their own. 

While often variable and unpredictable in brightness, LBVs 
are typically very luminous in the IR due to the significant 
amounts of dust formed in their outflows. Though hardly a 
homogeneous bunch, most LBVs are consistent with some 
sort of silicate dust in their spectra. Some, such as HR Car 
(lUmana et al.l 1200^ . show strong features from amorphous 
silicate grains, while others, like the Large Magellanic Cloud 
(LMC) LBV R71, exhibit features of crystalline silicates 
(iMorriset al.ll200l . Add itionally, a few LBVs like HDE 
31628 5 and HD 168625 (iMorris et all |200^ lUmana et all 
I2009h . show emission features from polycyclic aromatic hy¬ 
drocarbons (PAHs), though whether these PAHs are formed 
in the stellar outflow or are pre-existing in the i nterste llar 
medium (ISM) is unknown. In rj Carinae, ISmithI (1201 Ol) re¬ 
port that the IR spectra may be best fit with corundum dust 
(AI 2 O 3 ), which also shows a strong ^ 10 pm feature. We are 
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unable to identify any examples in the literature of IR spectra 
from CSM surrounding an LBV that unambiguously shows 
signs of graphite or amorphous carbon dust. 

Red supergiants (RSGs), typically very bright in the IR, 
also show signs of silicate dust in th e spectra of their C SM. 
In a sample of 19 RSGs in the LMC. IWoods et all (1201 11) re¬ 
port that amongst the 14 that show emission from warm dust, 
most show the 9.7 pm silicate feature, or other features of 
non-carbonaceous dust. They note that a few show weak PAH 
emission features, but that at least some of these may arise 
from the ISM. 

Other possibilities include any massive star that loses a sig¬ 
nificant amount of mass to stellar winds after evolving off 
the main-sequence. Yellow hypergiants are believed to be 
intermediate stages of stellar evoluti on with strong winds , 
between the RSG and LBV phases (iLagadec et all 1201 Ih . 
Again, though, the dusty shells seen around some of these 
stars, li ke IRC-h 10420, show str ong silicate features in their 
spectra (iHumphrevs et al.lll997^ . B[e] supergiants, massive 
15-7OM0), luminous, poorly understood early-type stars 
with slow stellar winds that form dusty circumstellar material 
may also be precursors to LBVs and/or WR stars. A Spitzer 
spectroscopic study of B[e] supergiants in the LMC shows 
that the dusty disks around these stars contain strong features 
of both crystalline and amorphous silicates, along with PAH 
features, but no carbonaceous dust (iKastner et al.ll200^ . 

An intriguing p oten tial progenitor is suggested by 
[Wesson et al.l (|2()10|) and iKochaneld (1201 111 based on archival 
Spitzer observations of the possible Type Iln SN 2008S. They 
rule out silicate grains and state that this absence of silicates 
may be evidence against a supergiant progenitor, favoring in¬ 
stead an “extreme” asymptotic giant branch (AGB) star. How¬ 
ever, this result should be interpreted with caution, as other 
authors report that SN 2008S is in fact a supernova impos - 
tor, and was more likely an LBV eruption (iSmith et alJl2009h . 
SN impostors and/or LBV eruptions may be relevant if, ul ti- 
mately, they explode as SNe Iln (see lMauerhan et £0120131) . 

Some authors have proposed a binary scenario for the pro¬ 
genitors of SNe Iln. Among binary systems, observations of 
carbon-rich stars AGB stars are consis tent with mass tr ans- 
fer from an evolved binary companion (lAoki et al.ll2002ll . In 
order to explain the apparent close t iming between the mass 
loss and SN events. IChevalied (l 20 1 2 h explored the possibility 
of mass loss driven by a common envelope involving a com¬ 
pact object and a massive star, finding that at least some SNe 
Iln could be explained this way. _ 

r] Car itself is a binary system, lining et aH (120051) report 
the likely detection of the companion star, rj Car B, and note 
that while a detailed stellar classification for this secondary is 
inconclusive, rj Car B could potentially be a Wolf-Rayet star. 
Neither the formation site nor the comp osition of the dust in 
rj Car is well-understood (iSmithll^OlOl) . Even the origin of 
the CSM is uncertain, as winds from either or both of the two 
stars could contribute. 

1 .2. SNe Iln in the IR 

Only a few SNe Iln have been observed and de tected at 
mid-IR (> 4.5 pm) wavelengths. iFox et alJ (120101) show 5- 
14 pm Spitzer IRAC and IRS observations of SN 2005ip, 
whose spectrum is clearly fit with a carbonaceous gr ain model 
with n o contribution from silicate dust required. Ivan DvkI 
( 1201 3h detected SN 1995N in archival data at 12 pm with 
the Wide-Field Infrared Survey Explorer, (WISE) and 24 pm 
with Spitzer, and showed that both silicate and graphite dust 


models have problems fitting the IR spectral energy distribu¬ 
tion (SED). Silicate grains were also ruled out for SN 2008S, 
which was detected at 24 pm by Spitze r, tho ugh its status as 
a supernova is questionable (see Section fLTI) . SN 2006jd was 
seen at 12 pm by the Wide-F ield Infrared Survey Explorer 
(WISE) (iStritzinger et al.ll2012^ without the silicate feature. 

SN 2010jl wa s observed from the ground in the 
optical/near-IR by iMaeda et alJ (l2013h . Their data is sensi¬ 
tive to the hottest component of the dust, and they find dust 
temperatures of > 1000 K. While they find equally good fits 
to this portion of the spectrum with carbonaceous and sili¬ 
cate grains, the temperatures required for the silicate grains 
are higher than the condensation temperature for those grains, 
and thus they favor carbonaceous grains. In this Letter, we 
report direct observations of the all-important mid-IR band 
using the Stratospheric Observatory for Far-Infrared Astron¬ 
omy (SOFIA), where we can detect the presence or lack of a 
silicate bump. We do not detect SN 2010jl, but the upper lim¬ 
its we derive rule out the presence of a silicate feature, and are 
w ell-fit by models of c arbon dust, confirming the conclusions 
of iMaeda et al.l (1201 3b . 

2. OBSERVATIONS 

SN 2010jl was discovered in UGC 5189A (distance 49 
Mpc) on 2010 Nov. 3.5 by iNewton & Puckett! (120101) . at a 
position of a = 09H2”53 .*337, <5 = +09 °29'42".13 (J2 ( )00.0) 
It v yas identified a s a Ty pe Iln by iBenetti et al.l (1201 ()l) 
and I Yamanaka et ^ (|2010|) . Multi-wavelength observations 
all point to a s trongly interacting SN with a dense CSM 
(iFransson et al.l 1201^ lOfek et aTri2010l) . Pre-SN Hubble 
Space Telescope images s how a luminous, blue point source 
at the position of the SN (ISmith et alJl201 ll) . but the data are 
insufficient to determine whether this source is a young clus¬ 
ter, a single luminous star, or a star caught during outburst. 
SN 2010jl shows significant evidence for late-time IR emis¬ 
sion. Some interpret th i s as evidence for newly formed dust 
(e.g. ISmith et al.1120121: iMaeda et al.ll2^ iGall et alJiMl 
iJencson et al.l 120151) . while others sugg est the presence of a 


pre-existing, unshocked dust shell (e. 

e.JAndrews et alJ 

20111 ; 

iFox et al.ll20131 IFransson et al.ll20l4 

Borish et al.ll2015 

). 


As part of an ongoing monitoring project of several Type 
Iln SNe at late times (PI; O. Fox), SN 2010jl has been ob¬ 
served in the mid-IR roughly 10 times with the Spitzer Space 
Telescope, all at the near-IR wavelengths of 3.6 and 4.5 pm 
(the SN occurred after the end of Spitzer’s cryogenic mis¬ 
sion). The SN is clearly detected in the Spitzer images, as 
we show in Figure [1] and has a flux of several mJy at both 
wavelengths. We also observed the SN with the FORCAST 
instrument on SOFIA using the 11.1 pm filter on 2014 May 
5-6 at an altitude of 38000 feet, with a total integration time 
of 6400 s. The images are taken in chop/nod mode with a 60" 
throw, and several hundred 30 s frames were coadded for a 
final data mosaic. The final mosaic has pixels 0.75" on a side 
with a resolution of approximately 2.7". 

Neither SN 2010jl nor the host galaxy were detected in our 
mosaicked image. The lack of detection of the host galaxy 
makes setting an upper limit on any point source easier, since 
we are background-limited and not confusion-limited. To ob¬ 
tain a 1 (T upper limit, we computed the standard deviation of 
the pixels in a 20 " x 20 " box centered on the source coordi¬ 
nates and converted this via the CALFCTR header parameter 
to get an RMS value in Jy/pix. We then multiplied this by 
the square root of the extraction area (4 pixel radius, 50 total 
pixels) centered on SN 2010jl. The upper limit obtained in 
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this method is 4.2 mJy, in good agreement with the 4.7 mJy 
value that we extrapolate from the SOFIA Observer’s Hand¬ 
book (Section 7.1.4) based on our exposure time. 

Since SNe are time-variable objects, we must compare our 
SOFIA observations with those most contemporaneous from 
Spitzer. We extract fluxes for four of the Spitzer observations, 
taken on 2013 Jan 30, 2013 Jun 29, 2013 Jul 3, and 2014 Jul 
9. We find no clear trend of systematic brightening or fading, 
so we average these four fluxes for each wavelength and use 
this average flux for our modeling, below. The average fluxes 
are 5.0 (1.4) and 5.69 (1.3) mJy for 3.6 and 4.5 p,m, respec¬ 
tively. The values in parentheses represent the standard devi¬ 
ation of the four measurements, which are significantly larger 
than Spitzer’s photometric uncertainties, and thus a very con¬ 
servative estimate of the uncertainties. 

3. DUST MODELING WITH THE NEW SOFIA UPPER LIMIT 

FI 1 modeled the Spitzer 3.6 & 4.5 pm emission from mul¬ 
tiple SNe Iln, considering a variety of origins and heating 
mechanisms for the dust. They rule out newly formed ejecta 
dust, as well as collisional heating of CSM dust by hot elec¬ 
trons behind the forward shock. They And that the most likely 
physical scenario is that a pre-existing dust shell is heated 
by UV and optical radiation from circumstellar interaction of 
the forward shock with dense material. The two short wave¬ 
length data points are insufficient to determine the type of dust 
present, and FI 1 considered both silicate and graphite grains. 
With the addition of the SOFIA upper limit at 11.1 /rm, we 
can compare models of silicate and graphite dust over a much 
broader wavelength range that includes the 9.7 pm silicate 
feature. 

The spectrum of a warm dust grain is simply the 
wavelength-dependent absorption cross-section of the grain 
multiplied by the Planck blackbody function. We con¬ 
sider several different types of grains. In addition to 
graphite and “astronomical” silicate grains, with absorption 
cross-sections derive d from bulk optical constants given in 
iDraine & Lee! d 1984ft. we also c onsider amorphous carbon 
grains (IRouleau & Martini 119911) . as well as other silicate 
grains, such as glassy and amorphous forms of both en statite 
and forsterite (iJd'ger et alJl2003t fPorschner et al.lll995h . For 
all grains, we assume a single grain size of 0.1 pm. This 
model is simplistic, but a more detailed model cannot be con¬ 
strained by only three data points, and the overall shape of 
the spectrum is not significantly cha nged by the addi t ion of 
more grain sizes. More importantly. iTemim & PwelB (1201 3h 
have recently shown that more physically realistic dust mod¬ 
els, which take into account a range of grain-size distributions, 
will lower the amount of mass inferred. Thus, our models 
represent an upper-limit to the amount of dust formed in the 
pre-existing CSM shell. 

We use a least-squares algorithm to fit the data for each 
grain composition. Since the flux values are of the same or¬ 
der of magnitude, we fit in linear space, and use only two 
free parameters: dust temperature and total mass. We ob¬ 
tain good fits to the SEP with models of both graphite (T 
= 553 K) and amorphous carbon (T = 616 K) grains, as we 
show in Figure |2] While we cannot further distinguish be¬ 
tween these two models, the difference is not particularly rel¬ 
evant to the goals of this paper. The only significant differ¬ 
ence between the graphite and amorphous carbon models is 
the amount of dust needed to fit the spectrum. Since amor¬ 
phous carbon has a higher emissivity than graphite over the 
range of ^ 1-20 pm, it requires less dust to achieve the same 


luminosity. Our graphite model requires 0.057 Mq of dust, 
while only 0.012 Mq of amorphous carbon is necessary (as¬ 
suming 49 Mpc distance). 

Most importantly, we And that no silicate grain model can 
adequately reproduce the data for all three points. We show 
the best-fit model for “astronomical silicates” in Figure |2l as 
fit to the two Spitzer points alone. This model requires a dust 
temperature of 863 K with a radiating dust mass of 0.024 Mq. 
The 11.1 pm model prediction from this model is roughly an 
order of magnitude higher than the upper limit of the flux den¬ 
sity from the SOFIA data. Even fitting a silicate model to the 
extreme uncertainty values of the Spitzer fluxes yields com¬ 
pletely nonphysical results. Other silicate models, such as 
“glassy” silicates, lead to even poorer fits. While we cannot 
rule out a small contribution from silicate grains (having only 
three data points makes this difficult to constrain), silicate dust 
does not dominate the SEP, and it is not required at all to ob¬ 
tain a good fit to the data. If the dust we see comes from a 
pre-existing shell of material, then this conclusion is at odds 
with the pre-SN dust observed in most proposed progenitor 
systems for SNe Iln. 

WR stars are one type of massive star t hat is known to pro - 
duce significant amounts of carbon dust (iWoods et alJl201 11 ) . 
However, since WR stars have blown off their hydrogen (and 
in some cases, helium) envelopes, it is generally accepted 
that WR stars would explod e as Type Ib/c SN e. An inter¬ 
esting suggestion is made by [Pwarkad^ (1201 Ih . who argues 
that at least some SNe Iln may explode after brief (< 10"^ yr) 
WR phases for their progenitors. This is quantitatively con¬ 
sistent with the model of Ell, who report that the progeni¬ 
tor system of SN 2006jd experienced a significant period of 
mass loss a few hund red years prior to explosion. Work by 
lAnderson et al.l (1201 2h found, by investigating the correlation 
of SNe Iln with the recent star formation history of the ex¬ 
plosion site, that SNe Iln trace recent, but not ongoing star 
formation, implying that most of these SNe do not arise from 
the most massive stars. 

There are other possibilites as well. The dust observed may 
be from newly-formed ejecta dust, and thus not indicative 
of the pre-SN outflow. A high optical depth or significantly 
larger grains th an are found in t he ISM might suppress the 
silicate feature (iLucv et al.lll991h . While we cannot say for 
certain that the lack of a silicate feature definitively means a 
carbon rich environment, the results here and in other SNe Iln 
are interesting and merit further study. 

4. SUMMARY 

The Type Iln SN 2010jl, easily seen in Spitzer 3.6 and 4.5 
pm, is undetected in 11.1 pm images from SOFIA. With the 
flux at the two Spitzer wavelengths and an upper limit at 11.1 
pm, we can characterize the dust composition. Only carbona¬ 
ceous grains (either graphite or amorphous carbon) can fit the 
spectrum; silicate grain models overpredict the 11.1 pm up¬ 
per limit by an order of magnitude or more. Our model fits 
suggest on the order of 0.01-0.05 Mq of dust radiating at a 
temperature of ^ 550-620 K. SN 2010jl joins a growing list 
of SNe Iln whose dust signatures do not show the 9.7 pm 
silicate feature. As El 1 argue, the dust in these SNe is con¬ 
tained in a pre-existing circumstellar shell, blown off by the 
progenitor system before the explosion. Thus, the dust that 
we observe can reveal something about the unknown nature 
of the progenitor. 

Most of the potential progenitor classes for SNe Iln are to¬ 
wards the massive end of the stellar spectrum: LB Vs, RSGs, 
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yellow hypergiants, and B[e] stars, all of which eject signif¬ 
icant quantities of material. However, they are all observed 
to produce silicate dust during their pre-SN lifetimes, not car¬ 
bonaceous. It is possible some stars may have a mixed chem¬ 
istry in their outflows, producing some amount of both silicate 
and carbon grains. This ha s been observed in some AGB stars, 
and IWilliams et al.] (1201 2|) inferred the presence of dust grains 
of both types in Spitzer observations of Kepler’s supernova 
remnant. Although Kepler is the remnant of a Type la, it is 
known to be interacting with a dense CSM from the progen¬ 
itor system, and could be a low mass analog of SNe Iln. An 
intriguing possibility is that of a binary system where one of 
the stars is a carbon-rich WR star. 

While we advise caution against drawing broad conclusions 
about the nature of SNe Iln from an indirect method and such 
a small sample size, we urge further study of these systems 
from the standpoint of both stellar evolution and dust nucle- 
ation in the atmospheres and outflows of massive stars. The 
James Webb Space Telescope will enable spectroscopic detec¬ 
tion of many SNe Iln in the mid-IR, allowing detailed studies 
of individual SN spectra, as well as the ability to do statistical 
studies on large samples of SNe Iln in the mid-IR. 

Facilities: SOFIA, Spitzer 
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Fig. 1. — Left: SOFIA 11.1 fj,m image of the field of SN 2010jl. Neither 
the SN nor the host galaxy are detected. Right: Spitzer 3.6 (blue) & 4.5 (red) 
image of the same region, clearly showing the SN, indicated by the tick 
marks, and host galaxy UGC 5189A. 
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Fig. 2.— Model fits to the data for graphite (solid black line), amorphous 
carbon (dashed line), and silicate dust (dot-dashed line). Spitzer 3.6 and 4.5 
pm fluxes are shown as squares, and the Icr SOFIA upper limit is marked with 
a triangle. Graphite and amorphous carbon dust models are fit to all three 
data points (assuming SOFIA upper limit is a detection), while the silicate 
dust model is fit only to the two Spitzer points. The best-fit silicate model 
overpredicts the 11.1 ^m emission by roughly an order of magnitude. 




